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The neurons in layer II of the medial entorhinal cortex are part of the grid cell network 
involved in the representation of space. Many of these neurons are likely to be 
stellate cells with specific oscillatory and firing properties important for their function. 
A fundamental understanding of the nonlinear basis of these oscillatory properties is 
critical for the development of theories of grid cell firing. In order to evaluate the 
behavior of stellate neurons, measurements of their quadratic responses were used to 
estimate a second order Volterra kernel. This paper uses an operator theory, termed 
quadratic sinusoidal analysis (QSA), which quantitatively determines that the quadratic 
response accounts for a major part of the nonlinearity observed at membrane potential 
levels characteristic of normal synaptic events. Practically, neurons were probed with 
multi-sinusoidal stimulations to determine a Hermitian operator that captures the quadratic 
function in the frequency domain. We have shown that the frequency content of the 
stimulation plays an important role in the characteristics of the nonlinear response, 
which can distort the linear response as well. Stimulations with enhanced low frequency 
amplitudes evoked a different nonlinear response than broadband profiles. The nonlinear 
analysis was also applied to spike frequencies and it was shown that the nonlinear 
response of subthreshold membrane potential at resonance frequencies near the 
threshold is similar to the nonlinear response of spike trains. 



Keywords: entorhinal cortex, stellate neurons, grid cells, quadratic sinusoidal analysis, frequency domain, 
nonlinear oscillations, resonance 



1. INTRODUCTION 

The stellate cells in layer II of the medial entorhinal cortex have 
long been noted for their oscillatory character (Erchova et al., 
2004) consisting of membrane potential oscillations (MPOs) and 
resonance properties (Engel et al., 2008; Giocomo and Hasselmo, 
2008; PastoU et al, 2012; Shay et al, 2012). More recently, it has 
been suggested that these cells participate in the grid-like fir- 
ing fields with regard to an animals position in space. In intact 
animals, grid cells have increased activity at particular locations 
representing a hexagonal grid. It is unclear how stellate cells 
along with pyramidal neurons participate in the grid cell net- 
works (Burgalossi and Brecht, 2014; Ray et al., 2014). It has been 
proposed that nonlinear stochastic current fluctuations from ion 
channels, rather than an internal periodic oscillator, is responsi- 
ble for this behavior (Erchova et al., 2004; Dodson et al., 2011). 
These neurons project to the hippocampus and are implicated in 
the activity of place cells that encode a single location. Thus, it is 
of some interest to understand in detail the nonlinear properties 
of stellate cells and how these are utilized in neural networks that 
compute spatial position. 

Quantitative information about the linear and quadratic 
behaviors can be obtained from a current clamp of the soma, 
which provides a basis to understand the just threshold behavior 



of these neurons in a particular neural network. We have used 
multi-sinusoidal stimulations to elicit stellate neuronal quadratic 
responses in the frequency domain, which are then encoded into 
a matrbc representing a Hermitian operator. The underlying the- 
ory, termed quadratic sinusoidal analysis (QSA), was published by 
Magnani and Moore (201 1 ) for voltage clamped neurons and in 
this paper, is applied to current clamp. The QSA not only quan- 
titatively characterizes experimental data from stellate cells, but 
it also provides an evaluation of the corresponding conductance 
based models. 

The method of frequency probing used here is based on a 
practical measurement technique, namely harmonic probing on 
Volterra kernels (Victor and Shapley, 1980; Boyd et al, 1983). 
Multi-sinusoidal stimulations enable nonlinear measurements of 
the neuronal response at harmonics 2/i and interactive frequen- 
cies \fi ± ^ I . These measurements constitute the coefficients of the 
QSA matrix, which represents an operator as an algebraic object 
similar to a Volterra kernel rather than just an array of numbers. 
It provides a precise signature of the nonlinear voltage dependent 
conductances and their particular representation on the dendritic 
and somatic membranes. 

Importantly, these studies show in detail how quadratic prop- 
erties of neurons are dependent on the frequency content of the 
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Stimulation and demonstrate that the nonhnear behavior accu- 
rately describes membrane potential responses at the level of 
normal synaptic activity (1-5 mV). 

2. MATERIALS AND METHODS 
2.1. TISSUE PREPARATION 

All experimental protocols were approved by the Boston 
University and University of Utah Institutional Animal Care and 
Use Committees. Horizontal sections of entorhinal cortex were 
prepared from 14- to 28 d-old LongEvans rats. All chemicals were 
obtained from Sigma (St. Louis, MO) unless otherwise noted. 
After anesthetization with isoflurane and euthanasia, brains were 
removed and immersed in 0°C artificial CSF (ACSF) consisting 
of the following (in mM): 125 NaCl, 25 NaHCOa, 25 D-glucose, 
2 KCl, 2 CaCl2, 1.25 NaH2P04, lMgCl2, and buffered to pH 
7.4 with 95/5% O2/CO2. Horizontal slices were cut to a thick- 
ness of 400 |xm (Vibratome 1000-I-; Vibratome, St. Louis, MO). 
Slices were incubated in a 32° C bubbled ACSF for 30 min before 
being cooled to room temperature (20°C). After the incubation 
period, slices were moved to the stage of an infrared, differential 
interference contrast-equipped microscope (Axioscope 2+; Zeiss, 
Oberkochen, Germany). All recordings were conducted between 



32 and 34° C. Solutions and preparation were identical to those 
described in Fernandez and White (2008). 

2.2. ELECTROPHYSIOLOGY 

Electrodes were drawn on a horizontal puller (P97; Sutter 
Instruments, Novato, CA) and filled with an intracellular solu- 
tion consisting of the following (in mM): 120 K-gluconate, 20 
KCl, 10 HEPES, 7 diTrisPhCr, 4 Na2ATP, 2 MgCb, 0.3 Tris-GTP, 
and 0.2 EGTA, buffered to pH 7.3 with KOH. Final electrode 
resistances were between 3 and 4M^2, with a range of mea- 
sured access values between 4 and 12M^2. All recordings were 
taken from the medial entorhinal cortex (MFC). Stellate cells were 
identified as neurons within layer II of the MFC exhibiting a 
large sag profile in response to hyperpolarizing current and hav- 
ing a peak between 2 and 7 Hz in their subthreshold impedance 
spectrum. Electrophysiological recordings were performed with 
a current-clamp amplifier MultiClamp 700A; Molecular Devices, 
Union City, CA), and data were acquired using custom soft- 
ware developed in MATLAB 2014 (MathWorks, Natick, MA) 
using the data acquisition toolbox. All additional details of the 
electrophysiological measurements are given in Fernandez et al. 
(2013). 



Stimulation Current (nA) 




-60 
-62 
-64 
-66 
-68 



0.4 0.8 
IVIembrane potential response at -64 mV 



time (s) 









s, ^ ^ ^1^^^ 

1 




1^ — . 

1 



0.4 0.8 time (s) 

Membrane potential response at -46 mV 




FIGURE 1 I Membrane potential response to multi-sinusoidal current 
stimulation with holding current adjusted to give —64 and — 46mV 
average membrane potentials. (A) Stimulation current with range of 



time (s) 



frequencies from 0.2 Hz to 2 kHz. (B) Membrane potential response 

at — 64mV. (C) Membrane potential response at —46 mV with spontaneous 

oscillations preceding the stimulation. 
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2.3. THEORETICAL SIMULATIONS 

The simulations were done by implementing Hodgkin-Huxley 
type conductance based models in MATHEMATICA 8 and 9 
(Wolfram Research, Champaign, IL, USA). The parameter esti- 
mation methods for all models is identical to that used for 
experiments on prepositus hypoglossi neurons in Idoux et al. 
(2008). The model structure consists of a soma and eight den- 
dritic compartments with uniform distributions for three voltage 
dependant ionic conductances : persistent sodium gNaP> potas- 
sium and hyperpolarization activated conductance ^h- In 
particular, since action potentials are not being simulated, the 
total sodium conductance is treated as one non-inactivating 
conductance gNaP- The gating variables are considered without 
power functions. Thus, the nonlinear behavior is essentially due 
to rate constants rather than power functions. Statistics were 
done with the PairedZTest of the Hypothesis Testing package of 
MATHEMATICA. 

2.4. QSA THEORY 

Linear systems are completely characterized by the linear super- 
position principle, which means that the response to a linear 
superposition of sine waves is a linear superposition of sine 
waves with the same frequencies but different amplitudes and 
phases. In contrast, the response of a nonlinear system can 
have frequencies not present in the stimulation. More precisely, 
if a multi-sinusoidal stimulation has frequencies f, and fj then 



the linear response will have frequencies and fj whereas the 
quadratic response will have additional harmonics 2/,, 2fj and 
interactive frequencies \fi ± ^- 1 . There also exist higher order inter- 
actions, such as fi+fj +fk, however the neurons studied in this 
paper mainly manifest quadratic nonlinearities. In particular, 
multi-sinusoidal stimulations require small amplitudes such that 
only linear and quadratic responses are significant, although the 
amplitudes must be large enough to overcome the background 
noise. 

A major obstruction to experimental measurements of non- 
linear responses is due to frequency overlaps. This means that 
two (or more) input frequencies can generate the same output 
frequency. For example, the input frequencies 1, 2, 3, 4 (in Hz) 
generate ambiguous output frequencies such as 3 — 1 = 4 — 2 (in 
Hz). A solution consists of choosing carefully the input frequen- 
cies without overlap up to the second order, higher orders being 
assumed negligible. An algorithm was written by Magnani and 
Moore (2011) to generate nonoverlapping frequencies. 

The Fourier transform applied to a multi-sinusoidal current 
stimulation /(f) generates Fourier coefficients at input frequen- 
cies^. The Fourier transform applied to the membrane potential 
V(t) generates Fourier coefficients at input frequencies/, as well 
as at harmonics 2/ and interactive frequencies \fi±fj\. In lin- 
ear analysis, the Fourier coefficients of lit) and V{t) can be 
used to compute the linear transfer function (impedance) Z[/jt] = 
V[fk]/I[fk]- Following Magnani and Moore (2011), we compute 




FIGURE 2 I Linear and quadratic analyses of a depolarized stellate 
neuron. (A) Zoomed measured membrane potential (solid), linear 
reconstruction (dashed) and linear + quadratic reconstruction (dotted). 
(B) Linear impedance at input frequencies (circle symbol at f,). (C) Quadratic 
output at harmonics (diamond symbol at 2f,), interactive sums (plus symbol 



at fj + fj) and interactive differences (inverted triangle symbol at |f/ — fj\]. The 
linear and quadratic analyses were truncated at 40 Hz in order to better 
illustrate the larger responses at the lower frequencies. Stimulation 
frequencies are 0.3, 1.3, 3.5, 6.4, 8.7, 12.1, 17.5, 23.7 30.8, 374, 42, 51.3, 
60.2, 74.9, 85.9 (in Hz). 
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FIGURE 3 I Effect of the membrane potential level on linear and 
quadratic responses of a stellate neuron at — 50mV (left column) 
and —69 mV (right column). The standard deviation (STD) is indicated in 
parenthesis. The frequencies are indexed as f(; where 

k er = {-N, -1 . +1 +W| and W = 15 denotes the number of 

stimulation frequencies. (A,B) Juxtaposed plots with respect to the input 
frequencies, namely the amplitudes of current stimulation /, linear impedance 
Z, and R function. A vertical bold line indicates the resonance frequency, 
which is f4 = 6.4 Hz for (A) and = 8.7 Hz for (B). The amplitude for each 
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plot ranges from zero to the indicated maximum. Each maximum coincides 
with zero for the next juxtaposed plot. (CD) Amplitudes of the coefficients of 
the QSA matrix Q, where the darkest rectangle represents the maximum 
value given by Max|Q| in the subtitle. Each matrix is indexed by r ordered 
from negative to positive numbers in the direction of the arrows. For 
example, column fe and row f_3 represents the frequency interaction 
fg - = 12.1 - 3.5 = 8.6 Hz. (E,F) Eigenvalues of the QSA matrix ordered by 
decreasing amplitudes. Stimulation frequencies are 0.3, 1.3, 3.5, 6.4, 8.7 
12.1, 175, 23.7 30.8, 374, 42, 51.3, 60.2, 74.9, 85.9 (in Hz). 



the quadratic transfer function 



Y'j- 



(1) 



where It is a time dependent vector encoding the multi-sinusoidal 
stimulation. 

The quadratic response can be reduced to a sum of squares 
through eigenanalysis of Q 



where ya = 1 and yij = 1/2 for i ^ j. The combinations fi + 
fj include all sums and differences between \fi\ and \fj\ when 
considering both positive and negative frequencies. The fre- 
quencies fk are indexed over the ordered set of integers F = 
{-N, ...,-1,-1-1, ...,-l-N}. 

Importantly, Magnani and Moore (2011) have shown that 
the complex matrix B can be turned into a Hermitian matrix 
Q, termed the QSA matrix, by row flipping of the coefficients 
Qi j = B^i j. Although a Hermitian matrix has complex coeffi- 
cients, its eigenvalues are real numbers, which are much easier 
to interpret physically as amplitudes in mV/nA^ . In this way, the 
quadratic part of the membrane potential can be expressed as an 
algebraic formula 

v2it)=i:Qit 



ier 

where d, are eigenvalues (mV/nA^) and |w,p are stimula- 
tion components (nA^) obtained from a linear transformation 
(Magnani and Moore, 201 1). In this way, the quadratic neuronal 
function can be interpreted as a set of quadratic filters where 
eigenvalues are amplitudes. When an eigenvalue is dominant and 
others are small, the quadratic neuronal function is approximately 
a simple square. 

Bode plots are useful to represent amplitude with respect to 
frequency. The QSA matrix can be represented in this way by 
summing the amplitude components of each column. However, 
unlike eigenanalysis, this induces information loss. The R sum- 
mation function (mV/nA^) is defined by 
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FIGURE 4 I Effect of the membrane potential level on linear and 
quadratic responses of stellate neuronal model simulations at —49 mV 
(left column) and — 57mV (right column). (A,B) Juxtaposed plots with 
respect to the input frequencies, namely the amplitudes of current 
stimulation /, linear impedance Z, and R function. The resonance frequencies 
are respectively fs = 8.7 Hz (A) and fg = 12.1 Hz (B). (C,D) Amplitudes of the 
coefficients of the QSA matrix Q with the maximum value given by Max|Q|. 
(E,F) Eigenvalues of the QSA matrix. The model parameter values (in units 
based on mV and nA) are as follows : membrane capacitance 



Csoma = 0.0000542, maximal conductances geak = 0.0005, gK = 0.023, 
gi^gP = 0.0024 and gn = 0.014; the reversal potentials V^eak = ^55, 
Vk = -87, VNap = 77 and Vh = -43. The functions an, Pn, am, Pm, a/j, Ph 
depend on the variable V and their description is given by Idoux et al. (2008) : 
Vn = -35, Sn = 0.045, tn = 0.7543, Vm = -38, Sm = 0.06, tm = 0.000150, 
Vf, = —51.65, Sf, = —0.06 and t/, = 0.05. Finally, the electronic length is 
eiength = 0.50 and the ratio of the dendritic area to the soma area is 
A,stio = 4.5. Stimulation frequencies are 0.3, 0.4, 1.3, 1.5, 3.5, 4.0, 6.4, 8.7 
12.1, 175, 21.6, 23.7 30.8, 374, 42.0, 51.3, 56.9, 60.2, 74.9, 85.9 Hz (in Hz). 



J^j = ^IQ.;l (2) 
ier 

Alternatively, a modified function R' with units homogenous to a 
linear transfer function (mV/nA) is defined by 

Intuitively, the R and R' functions evaluate the sum of the ampli- 
tudes at all frequencies for which they interact. 

Figure 1 illustrates the membrane potential response of a 
stellate neuron to a multi-sinusoidal current stimulation at 
two membrane potential levels. Membrane potential oscilla- 
tions (MPOs) are clearly observed before and after the stim- 
ulation at the depolarized level. Also, the response shows 
larger oscillations during the depolarization. This neuron has 
a strong resonance in the same frequency range as the MPO's 
(Schreiber et al, 2004; Engel et al, 2008; Yoshida et al, 
2011). 



Figure 2 illustrates the linear and quadratic analyses of a 
depolarized stellate neuron. Figure 2A compares the measured 
membrane potential in the time domain to the linear recon- 
struction and linear -|- quadratic reconstruction. Clearly, for 
these signal levels, linear analysis is insufficient and quadratic 
analysis is required to capture the nonlinear neuronal behav- 
ior. Figures 2B,C show the linear impedance (mV/nA) and the 
quadratic output (mV) in the frequency domain. The Fourier 
components of the linear impedance occur at the stimulation 
frequencies, whereas those of the quadratic output occur at sec- 
ond order harmonics and interactive frequencies. Clearly, the 
quadratic output shows a strong overlap with the linear resonance 
frequencies. 

The QSA theory applied to neurons does not rely on any 
assumption about their morphology or ion channel distribution. 
In particular, the complexity of Hodgkin-Huxley based multi- 
compartmental models with their computational overhead is 
avoided. QSA analyses were done with MATLAB on both model 
simulations and experimental data, allowing precise comparisons 
of the quadratic responses. 
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3. RESULTS 

3.1. EFFECT OF MEMBRANE POTENTIAL 

Figure 3 illustrates the linear and quadratic responses of a stel- 
late neuron for two membrane potentials —50 and —69 mV 
(left and right columns respectively). The stimulation frequen- 
cies were identical to those of Figure 2 and are indexed as ft 
where fc e F = {—15, . . . , —1, -1-1, . . . , -1-15}. The upper panels 
Figures 3A,B represent a juxtaposition of the amplitudes for the 
stimulations I, linear impedances Z, and -R functions with respect 
to the input frequencies. The -R functions are plotted as Bode 
plots in the same way as the impedances although the ordinate 
units are different. The maximum of each R function is close 
to the impedance resonance frequency. Statistics were calculated 
for a group of six stellate neurons. The maximum amplitude of 
the QSA matrix increased from 275 to 715 mV/nA^ (p = 0.0004) 
for a membrane potential change of -|-7mV in the range —65 
to -48 mV. 

As explained in previous publications (Magnani and Moore, 
2011; Magnani et al., 2013), the QSA matrix provides a com- 
plete description of the quadratic response as ratios between 
output and input coefficients (Equation 1). In contrast, the coded 
points in Figure 2C show the quadratic measurements of the 
output without showing the frequency interactions ifi,fj). The 
QSA matrix was constructed from the Fourier coefficients of the 
data, then decomposed by eigenanalysis to compute its eigen- 
values. Each cell of the QSA matrix represents the amplitude of 
the voltage response at two interactive frequencies divided by the 
amplitude of the current at these two frequencies (Equation 1). 
The matrix plots have many symmetries, which reflect particu- 
lar algebraic properties of the underlying neural operator. The 
matrix plot is indexed by the set F ordered from negative to 
positive numbers in the direction of the arrows, that is to say 
/_15, . . . . . . ,/-|-i5. Each cell at abscissa /J and ordinate 

fj encodes the ratio between the membrane potential at f, + fj and 
the current at // and fj. In particular, the white diagonal encodes 
interactions fi + = 0, thus the DC is set to zero. The other 
diagonal/; + fi encodes the harmonics 2/,. In Figures 3C,D, the 
QSA amplitudes appear concentrated at lower frequencies near 
the center of the matrix, which is consistent with the maximum 
of the R function. 

The impedance and resonance frequency are much less depen- 
dent on the membrane potential than the QSA matrbc and R 
function. This suggests that the quadratic neuronal function espe- 
cially encodes nonlinear voltage dependent ionic conductances. 
The effect of the membrane potential on stellate neurons is pro- 
nounced for all nonlinearities, namely the amplitudes of the QSA 
coefficients, the eigenvalues and the _R functions. 

In this and all subsequent figures, the -R function (mV/nA^) is 
juxtaposed on the linear impedance (mV/nA), which in turn is 
juxtaposed on the stimulation amplitude Fourier spectrum (nA). 
Although the R function is a non reversible reduction of the QSA 
matrix, it provides a practical way to compare the linear and 
quadratic behaviors at input frequencies. It can be observed, in 
Figures 3A,B, that the R function has a resonance frequency range 
comparable (but not identical) to the linear case. 

At —50 mV (Figure 3, left column), the QSA matrix gives more 
detail on frequency interactions showing enhanced amplitudes in 



the centered square delimited by |/±4| = 6.4 Hz, namely those 
that involve the resonance frequencies. At the more hyperpo- 
larized membrane potential —69 mV (Figure 3, right column), 
the QSA matrix shows lower amplitudes. However, there is a 
peak for the harmonics of the lowest input frequency 2/i = 
0.6 Hz. Moreover, the quadratic response is enhanced around the 
centered square delimited by the slightly higher resonance fre- 
quency/5 = 8.7 Hz. Thus, the shift in frequency of the nonlinear 
responses with membrane potential level is similar to the voltage 
dependence of the linear resonance frequencies (Shay et al., 2012). 

The eigenanalysis of the QSA matrix reveals that the non- 
linear function is concentrated in a single dominant eigenvalue. 
This suggests that the neuronal processing consists of a single 
nonlinear-linear unit as opposed to a parallel combination of sev- 
eral units (see Magnani and Moore, 2011). Dominant eigenvalues 
are frequently observed, however there are generally multiple sig- 
nificant eigenvalues. Figure 4 illustrates that a simplified model 
of stellate neurons is able to capture the nonlinear behavior 
described by experiments of Figure 3. 

3.2. LINEAR DISTORTIONS 

Previous experiments by Haas and White (2002) suggested that 
the Fourier responses of stellate neurons to multiple frequencies 
were dependent on the frequency content of the stimulation, and 
furthermore these responses were different from that obtained 
when using single sine waves. We have explored this issue in 
more detail by comparing the responses of stellate neurons to 
stimulations containing many frequencies vs. QSA stimulations 
containing nonoverlapping frequencies. The presumed linear 
responses in the mV range to many frequencies show marked 
distortions compared to almost no distortions using nonoverlap- 
ping frequencies. The term, distortions, is used here to describe 
an irregular function with larger variations than generally seen in 
our linear impedance functions. 

Figures 5A,C,E (left column) illustrates distortions in the 
frequency responses from stellate neurons due to nonlinear 
responses at overlapping frequencies. The upper trace (A) shows a 
marked distortion at —53 mV of linear responses stimulated with 
a very large number of frequencies near threshold. The middle 
trace (C) shows a significant reduction in the distortion of the lin- 
ear response for a more hyperpolarized membrane potential with 
the same number of overlapping frequencies. Finally, the bot- 
tom trace (E) shows an almost undistorted linear response using 
nonoverlapping frequencies (up to the second order) despite a 
greater membrane potential STD. 

Clearly, stimulations with overlapping frequencies generate 
interactive frequencies that overlap the input frequencies and 
thereby cause a distortion in the linear response. These effects nat- 
urally increase with the number and amplitudes of the particular 
frequencies. Similarly, Figures 5B,D,F (right column) illustrates 
similar distortions from a neuronal model. The upper panel (B), 
near threshold, shows that the linear response is significantly dis- 
torted due to the presence of nonlinear interactive frequencies 
for potential excursions in the mV range. The middle panel (D) 
shows less but distinct distortions at a lower STD and a more 
hyperpolarized membrane potential, which is consistent with 
the experimental results (left middle panel). The bottom panel 
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FIGURE 5 I Distortion of linear impedance (mV/nA) by tlie stimulation 
profiles for experimental data and model simulations. Tine average 
membrane potential is given along with STD values of membrane potential. 
The standard deviations are given in parenthesis. (A,C) Experimental data : the 



200 



Model Impedance (mV/nA) 
Vm=-53mV(±0.6) 




200 



5 10 
Vm=-56mV(±0.26) 




5 10 
Vm=-55mV(+0.06) 



F 

200 
150 
100 
50 



same overlapping frequencies in 0.3 — 1000 Hz with 0.1 Hz intervals. 
(E) Experimental data : the nonoveriapping frequencies are the same as 
Figure 4. (B,D,F) Model simulations : the same overlapping frequencies in 
0.3 — 500 Hz with 0.1 Hz intervals. The model parameters are given in Figure 7. 




(F) shows an almost undistorted linear response with the same 
overlapping frequencies for a still lower STD membrane potential. 

Generally, if the output amplitude STD using many overlap- 
ping frequencies is less than 0.5 mV, then the linear responses 
have minimal distortion. Undistorted linear responses can be 
obtained with fewer nonoveriapping frequencies for much higher 
STD membrane potential responses, typically up to 4 mV. Thus, 
a larger number of overlapping frequencies evoking responses in 
the mV range induce significantly more nonlinear effects than 
fewer nonoveriapping frequencies. Nonstationary signals, such as 
the chirp or the zap, are not optimal for linear analysis since their 
Fourier components distribute over all frequencies and dramat- 
ically generate frequency overlaps. Nevertheless, linear analysis 
is still possible for small stimulation amplitudes (Erchova et al., 
2004; Schreiber et al, 2004). 

3.3. LOW PASS STIMULATION FILTERING 

The above results suggest that the effect of the frequency con- 
tent of the stimulation on the true nonlinear response should 
be distinguished from the linear distortions. Experiments were 
done with a low pass filtered stimulation using nonoveriapping 
frequencies in which the high frequency amplitudes are reduced. 
Nonoveriapping frequencies assure that the linear responses 
are not distorted in order to accurately measure the nonlinear 
behavior. Since stellate neurons have resonance frequencies 



around 10 Hz or lower, it would not be surprising that these fre- 
quencies compared to high frequencies would have the largest 
amplitude responses and thus can easily evoke the nonlinear 
behavior. 

The nonlinear responses induced by differently filtered 
Gaussian white noise inputs are likely to be stimulation depen- 
dent. This not only leads to nonlinear frequency interactions 
dependent on the frequency content of the stimulation, but linear 
analysis can also be altered since Fourier components at individ- 
ual input frequencies would be contaminated by nonlinearities. 
Surprisingly these effects can occur at relatively small amplitudes 
of membrane potential responses, namely in the mV range. Thus, 
empirically, responses of neurons at their stimulating frequencies 
can be dependent on the frequency content of the stimulation, 
although a true linear response is not. 

Figure 6 shows that the nonlinear responses of a stellate neu- 
ron to various low pass filtered stimulations are remarkably 
different despite identical linear behaviors. The cutoff frequen- 
cies are not the same between left and right columns, showing 
that the increase of high frequency interactions are related to high 
frequency filtering. Statistics were calculated for a group of three 
stellate neurons comparing the effect of the two low pass stimula- 
tion filters. The maximum amplitude of the QSA matrix increased 
from 1229 to 4204 mV/nA^ (p = 0.019) for a membrane potential 
change of —60 to —55 mV. 
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FIGURE 6 I Effect of stimulation profile on stellate neurons. Left column 
(A,C,E) Weakly low pass filtered stimulation / in dotted lines. Right column 
(B,D,F) Strongly low pass filtered stimulation / in dotted lines. The resonance 
of the impedance (Z) is fg = 4.0 Hz for both stimulations (A,B). The 



attenuation of high frequencies in the stimulation of the second column 
induces enhanced quadratic responses as indicated by the R functions in the 
top row (A,B) as well as the QSA matrix (CD) and eigenvalues (E,F). The 
stimulation frequencies are the same as in Figure 4. 



The coefficients of the QSA matrix are defined as the ratios 
between an interactive output and the product of the correspond- 
ing inputs for each pair of input frequencies. Thus, although the 
amphtudes of the high frequency responses decrease when the 
stimulation is filtered, the above ratio tends to increase. Therefore, 
the quadratic function appears sensitive to the frequency content 
of the stimulation even for membrane potential responses in the 
mV range. 

Increases in the R functions at high frequencies are clearly 
observed, however the QSA matrix gives more detail about the 
interactive frequencies. The QSA matrix of the left column shows 
peak amplitudes for the harmonics of input frequencies near 
the linear resonance frequency, such as 2/5 = 8 Hz. Moreover, 
there are significant interactive differences between high input 
frequencies such as/19 — /18 = 14.7 Hz. The QSA matrix of the 
right column shows greatly enhanced high frequency interac- 
tions in response to the high frequency filtered stimulation. 
Although the lower frequency interactions near the resonance 
are not apparent, their amplitudes are similar as can be seen 
from the R functions. The complexity of these responses is also 
reflected through eigenanalysis that turns a single dominant 
eigenvalue to multiple eigenvalues when the stimulation is fil- 
tered. Figure 7 illustrates similar results for model simulations. 
It will be shown below that the increases of the responses at 
interactions between high frequencies are quite sensitive to the 



relative amplitude ratio of the low vs. high frequency stimulation 
content. 

The stimulation dependence of the nonlinear response is a 
good example of the dramatic difference between linear and non- 
linear behaviors even at small signal amplitudes. Linear responses 
are independent of stimulation amplitude for all frequencies. 
However, we have shown that interactive differences between high 
frequencies are not proportional to the stimulation amplitude 
when sufficiently low pass filtered. The quadratic transfer func- 
tion is defined as the ratio between an interactive output and 
the product of two inputs. When the enhanced interaction cor- 
responds to a difference between high frequencies, the ratio tends 
to be maximized because the numerator reflects a significant non- 
linear response and the denominator is small due to the low pass 
filtering. 

3.4. BAND PASS STIMULATION FILTERING 

The above results demonstrate that the frequency content of the 
stimulation plays an important role in the nature of the nonlinear 
response, but generally has minimal effects on the linear behavior 
if the membrane potential responses are less than 0.5 mV. In order 
to determine if nonlinear responses have limiting small signal 
responses that are independent of the content of the stimulation, 
model simulations were done by applying a Gaussian-like win- 
dow to the stimulation in the frequency domain. In this way, the 
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FIGURE 7 I Effect of stimulation profile on stellate neuronal model 
simulations. Left column (A,C,E) Weakly low pass filtered stimulation / in 
dotted lines. Right column (B,D,F) Strongly low pass filtered stimulation / 
in dotted lines. The resonance of the impedance (Z) is fe = 4.0 Hz for both 
stimulations (A,B). The attenuation of high frequencies in the stimulation 
of the second column induces enhanced quadratic responses as indicated 
by the R functions in the top row (A,B) as well as the QSA matrix 
(CD) and eigenvalues (E,F). The stimulation frequencies are the same as 
in Figure 4. The parameter values (in units based on mV and nA) are as 



follows : membrane capacitance Csc 



-. 0.0000542, maximal 



conductances gieak = 0.0005, = 0.0627, gNaP = 0.0065 and 
gH = 0.002926; the reversal potentials Vjeak = -55, Vk = -87, V^^p = 77 
and Vh = —43. The functions «„, fin, am, fim, a/i, Ph depend on the 
variable V and their description is similar to those published by Idoux et al. 
(2008) : v„ = -38, = 0.035, t„ = 0.75, Vm = -38.76, Sm = 0.046, 
fm = 0.000150, i/h = -51.65, = -0.01 and f/, = 0.5567. Finally, the 
electronic length is eiength = 0-50 and the ratio of the dendritic area to the 
soma area is Aatio = 4.5. 



Stimulation shows a peak amplitude at a specific frequency (like a 
resonance). 

Figure 8 illustrates that, in contrast to the linear impedance, 
if the stimulation peak is at a low frequency (panel A, near 
the resonance), then the maximum of the QSA matrix is much 
greater than for a stimulation peak at a higher frequency (lower 
row, beyond the resonance). Both stimulations evoke potential 
responses in the mV range and have identical linear impedances. 

Discrete peaks can also be observed for the QSA matrix in 
Figures 8B,D. The panel B shows discrete peaks close to the har- 
monics of/4 = l.SHzas well as high frequency interactions in the 
corners. The panel D also shows discrete peaks clustered around 
harmonic frequencies of fj = 6.4 Hz and /g = 8.7 Hz. In sum- 
mary, Gaussian-like window applied near the resonance (panel A) 
evokes output frequencies before and after /g = 8.7 Hz. Gaussian- 
like window applied after the resonance (panel C) leads to output 
frequencies that are near the resonance. Therefore, nonlinear 
responses involve a different frequency range compared to the 
linear resonance frequency. 

Figure 9 illustrates that the quadratic response to very small 
stimulations is almost independent of the frequency content. In 



particular, both linear and quadratic functions are nearly identi- 
cal for flat stimulation vs. low peak stimulation, as well as higher 
peak stimulation (not shown). Thus, the enhanced corners for 
low peak stimulation almost disappear if the membrane poten- 
tial responses have STD values of 0.1 mV or less. In particular, the 
QSA matrix has a remarkable sensitivity to stimulations with low 
frequency content, especially near the resonance of the neuron. 
Indeed, Figure 9A shows a larger R function for such nonuniform 
stimulation profile. 

From a mathematical point of view, the limit value of the QSA 
matrbc for very small stimulations is consistent with local anal- 
ysis and represents a quadratic transfer function as a Volterra 
kernel. For larger stimulations, the quadratic response becomes 
dependent on the stimulation and the ratio (Equation 1) can- 
not be interpreted as a transfer function. However, the non- 
linear output is still a quadratic signal without higher order 
frequency contamination (Magnani and Moore, 2011). In fact, 
larger low peak stimulations evoke the low frequency range near 
the resonance, in which interactions l/j- — fj \ corresponding to 
the enhanced corners are generated. Such an interpretation is 
supported by the enhanced output amplitudes of the interactive 
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FIGURE 8 I Nonlinear effects of a Gaussian-like window applied to the 
stimulation Fourier components. Upper row (A,B) Gaussian-like window 
applied to the stimulation Fourier components near the resonance 
fs = 8.7 Hz. This induces enhanced frequency interactions between input 



frequencies near the resonance. Lower row (C,D) Gaussian-like window 
applied to the stimulation Fourier components three fold above the 
resonance fs = 8.7 Hz. The QSA matrix appears much smaller than for a 
stimulation peak near the resonance. 



difference frequencies in Figure 2C (inverted triangles). This may 
increase the quadratic response at these interactions, relative to 
the stimulation. 

These nonlinear effects are present at the usual membrane 
potential excursions in the mV range and appear amplified 
for stimulations tuned to the preferred frequency response 
of a particular neuron. Interestingly, other constant stimu- 
lation profiles (not represented) lead to quadratic functions 
independent of stimulation amplitudes for membrane poten- 
tial responses up to a few mV. For stimulation profiles with 
enhanced low frequencies, this independence of stimulation 
amplitude occurs only for a tenfold lower response ampli- 
tude (10-100 microvolts), which again emphasizes the sensitiv- 
ity of the nonlinear responses to the frequency content of the 
stimulation. 

These results are consistent with experimental data of Figure 6 
showing that low frequency stimulation leads to an enhance- 
ment of high frequency interactions for both QSA matrix and 
R function. Figure 8 indicates that high frequency interactions 
elicited by low frequency stimulation can nearly disappear if the 
stimulation peak is shifted to higher intermediate frequencies. In 
both cases, the linear responses remains identical despite major 
differences in the nonlinear behavior. 



3.5. DENDRITIC STIMULATION 

Since cable properties of dendrites induce a signal filtering, 
the marked effect of the stimulation content on the nonlinear 
response could play a role in how dendrites process signal inputs. 
It can be expected that distal vs. proximal synaptic inputs to single 
neurons would lead to different responses in the soma. Figure 10 
illustrates such effects with model simulations. It compares the 
soma potential response to a current stimulation in the soma (left 
column) vs. a stimulation of the most distal dendritic compart- 
ment (right column). In both cases, the stimulation profile is 
flat but the represented current I is not the stimulation but the 
computed current flowing across the soma membrane. 

In order to compare the two approaches, a somatic impedance 
was calculated by determining the ratio of the soma potential and 
the current flowing across the soma membrane, rather than the 
total injected current typically used for either point or transfer 
impedances. This procedure determines just the impedance of the 
somatic compartment in isolation, which is nearly independent of 
the stimulation profile, as shown in Figure 10. At smaller stimula- 
tions, the two impedances converge to identical functions. Thus, 
this approach can be used to analyze the effect of dendritic fil- 
tering of the injected current from the dendrite on the nonlinear 
responses of just the soma. 
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FIGURE 9 I Equivalence of nonlinear responses for different stimulation amplitude stimulation profile. The linear and nonlinear responses are nearly 
profiles at very small amplitudes. Upper row (A,B) Gaussian-lil<e window the same, however the maximum of the R function is about two fold higher 
applied to the stimulation Fourier components. Lower row (C,D) constant for the Gaussian-like stimulation. 



In general, the computed impedance for the somatic com- 
partment is higher than the impedance of the soma plus the 
dendrite. This occurs because the impedance of a single somatic 
compartment is greater than the soma with attached dendritic 
compartments (not shown). More precisely, a significant por- 
tion of the injected current passes down the dendrite leaving a 
smaller current for the soma, which has a minimum value near 
the resonance frequency, at which the impedance is by definition 
maximal. Thus the smaller somatic current with the same soma 
potential gives a greater impedance. The total injected current 
would be the sum of the somatic current and remaining current 
flowing down the dendrite. 

When a flat stimulation is injected at the end dendritic com- 
partment (Figure lOB), the somatic current has a minimum near 
the resonance frequency and the somatic impedance is nearly the 
same to that determined with somatic stimulation. The simula- 
tions show that the nonlinear responses (QSA matrix, R function) 
are greater for the distal stimulation (right column) due to the 
dendritic filtering, which can be seen to be a Gaussian-like peak 
around 15 Hz. 

When stimulating the end dendritic compartment, the other 
compartments significantly reduce the high frequency content of 
the current reaching the soma (Figure 10, right column). This 
current increases from a minimum at the resonant frequency to 



a maximal value and then progressively decreases with frequency 
as would be expected from dendritic filtering. These effects are 
the consequence of the progressive filtering of the membrane 
potential by each dendritic compartment toward the soma. Each 
individual dendritic compartment shows similar scaled linear 
behavior, which is dependent on the potential in each compart- 
ment as determined by the cable structure including the voltage 
dependent channels. 

The QSA matrix for the dendritic stimulation shows enhanced 
corners that are similar to those observed for filtered stimulations 
applied to the soma in Figures 6, 8. However, the dominating 
responses at the intermediate harmonics and interactive differ- 
ence frequencies are comparable for both somatic and dendritic 
stimulations. Decreasing the stimulation amplitudes leads to the 
same limit functions for impedances and QSA matrices (not 
shown) as found for the Gaussian-like stimulations of Figure 9. 

Figure lOD also shows peaks for the QSA matrix at harmon- 
ics and high frequency interactions (around l/e = 8.0 Hz and 
/20 — /i9 = 11 Hz). These are similar to those observed for the 
Gaussian-like window of Figure 8, which has peaks clustered 
around harmonic frequencies of/7 = 6.4 Hz and/s = 8.7 Hz. 

Thus, the dendritic filtering is simulated by a Gaussian-like 
stimulation applied to the soma, which is related to the effect of 
a filtered stimulation applied to the soma as in Figures 6, 7. In 
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FIGURE 10 I Effect of distal vs. proximal flat stimulation on somatic 
nonlinearities. The linear and quadratic functions were calculated using 
only the current passing across the soma membrane, as opposed to the 
total injected current. Thus, the analyses apply to the soma independent 
of the location of the stimulating electrode. Left column 
(A,C,E) Responses to the flat stimulation of the somatic compartment. 
Right column (B,D,F) Responses to the flat stimulation applied at the 
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most distal compartment. (A,B) Juxtaposed plots with respect to the 

input frequencies, namely the amplitudes of current stimulation /, linear 

impedance Z, and R function. The resonance of the impedance is 

fg = 4.0 Hz for both stimulations. (CD) Amplitudes of the coefficients of 

the QSA matrix Q with the maximum value given by Max|Q|. 

(E,F) Eigenvalues of the QSA matrix. All parameters of the model neuron 

are the same as in Figure 7. 



summary, these results support the fact that the enhancement of 
low frequency amplitudes increases nonlinear responses. 

Increasing the electronic length leads to somatic current pro- 
files similar to the filtered stimulations of Figures 6, 7 and QSA 
matrices that are quite similar. From a physiological point of view, 
these nonlinear effects due to filtering are very much dependent 
on the precise dendritic location of the stimulation and appear to 
be present at all input amplitudes. However, these nonlinearities 
are not significant until the soma potential responses are in the 
mV range where they also contribute to a distortion of the linear 
behavior if the frequency content of the stimulation is overlap- 
ping. Clearly, higher order nonlinearities are also generated with 
very large stimulations. 

3.6. SPIKE FREQUENCY MODULATION 

The subthreshold nonlinear responses are markedly dependent 
of the membrane potential and essentially reach their maximum 
values just below the threshold. This would suggest that the 
nonlinear behavior at just below the threshold membrane poten- 
tial could provide a reasonable estimation of the suprathresh- 
old action potential response. In order to directly determine 
suprathreshold responses, a QSA analysis was done on action 



potential responses using three or four stimulation frequencies 
near the resonance, namely 0.3, 1.3, and 3.5, or 4, 6.4, 8.7, and 
12.1 Hz, the latter encompassing the resonance range of frequen- 
cies. Both linear and quadratic analyses were done using a Fourier 
analysis of unit spike events constructed from their peak values, 
which extracts the expected low frequency responses. 

Figure llA shows modulated action potential responses to a 
multi-sinsuoidal stimulation containing three frequencies (0.3, 
1.3, and 3.5 Hz) that are all below the resonance. The under- 
lying traces show the linear and quadratic membrane potential 
reconstructions for just these frequencies, which are sufficiently 
low to minimize the contamination of the Fourier transform 
by the shape of action potentials. Figure 1 IB represents linear 
and quadratic analyses of the membrane potential for the same 
neuron stimulated just below threshold. Since there are few fre- 
quencies, the superimposed reconstructed linear and quadratic 
responses have slight differences (Figure 11 A, dashed linear vs. 
dotted quadratic lines). However, the quadratic responses remain 
significant (Magnani and Moore, 2011) to encode fundamen- 
tal nonlinear properties of the neuron that cannot be predicted 
by linear analysis. The linear impedance shows three rising val- 
ues that precede the usual resonance of stellate neurons. The 
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FIGURE 11 I Linear and quadratic membrane potential responses for 
pre-resonance stimulation frequencies (fi = 0.3, f2 = 1.3 and fi = 3.5 Hz). 

(A) IVIembrane potential responses (solid line) with modulated action potentia 
responses. The dashed and dotted curves show the linear and quadratic 



r={-N,...,-i,+i,...+N) 



responses of the underlying membrane potential for the three stimulating 
frequencies. (B) Linear and quadratic analyses for a subthreshold response 
displaced by 0.5 mV without action potentials. The corresponding reconstructed 
time domain analyses are almost identical to those shown in the first row. 



QSA matrix and the R function indicate that both low and 
higher frequencies near the harmonics have the greater inter- 
active responses, which are also chararactized by a dominant 
eigenvalue. 

The time instants of each action potential of Figure llA were 
converted to unit impulses (UI), as illustrated in Figure 12A. The 
dashed and dotted traces show the linear and quadratic analy- 
ses computed from the Fourier transform of the unit impulses. 
Figure 12A shows the waveforms of both the linear and quadratic 
analyses of the unit impulses, which are similar to those of the 
membrane potential of Figure 1 lA. The QSA matrix and the R 
functions are slightly different, showing a greater amplitude for 
the highest frequency interactions, mainly the highest harmonic 
frequencies 2/3 = 7 Hz. It is evident from Figure 12A that, when 
the frequency of spiking modulation is high, namely the dotted 
quadratic response has a much greater amplitude than the dashed 
linear line. Also, both the membrane potential and unit impulses 
analyses have a dominant eigenvalue, which suggests that they 
have features in common for their nonlinear behavior. 

Since the Fourier analysis of a spike train contains many high 
frequency components unrelated to the lower modulated spiking 
frequencies, it is useful to construct an instantaneous frequency 
(IF) curve from the spike times. A slightly Gaussian filtered 
IF curve is shown in Figure 12C as a solid trace for the same 
time range illustrated for the unit impulses in Figure 12A. The 
dashed linear and dotted quadratic traces are the reconstructed 
responses of the solid trace. The linear and quadratic analyses 



are shown in Figure 12D and are similar to the analysis of the 
unit impulses in Figure 12B except for the details of the QSA 
matrix. The QSA matrix of unit impulses shows high frequency 
interactions at sums + fj and harmonics 2/,. The QSA matrix 
of instantaneous frequency curve has increased amplitudes of 
high frequency interactions at differences \fi — fj\- These various 
effects are also reflected in the time domain of the reconstructed 
responses (first row compared to third). 

Figure 13 illustrates linear and quadratic analyses of the mem- 
brane potential (Figure 13A) measured without action potentials 
compared to unit impulses (Figure 13B), for a stimulation based 
on four frequencies (/l = 4,/2 = 6.4, /s = 8.7 and/4 = 12.1 Hz) 
that encompass the resonance of a stellate neuron. It is apparent 
that both the linear and nonlinear behaviors are quite compara- 
ble in this frequency range. The quadratic responses associated 
with the middle two frequencies are enhanced for both the mem- 
brane potential and unit impulses analyses. In general, the QSA 
matrix of unit impulses shows a more restricted set of interme- 
diate frequency interactions than observed for the QSA matrix 
of the membrane potential, as can be observed on Figure 13. 
Nevertheless, the similarity of the linear and quadratic analy- 
ses of the subthreshold membrane potential and unit impulses 
are striking. This suggests that the suprathreshold modulated 
spiking behavior is reasonably well approximated by the sub- 
threshold membrane potential nonlinearity, just below threshold. 
As threshold is approached, the subthreshold responses take on 
a more nonlinear character, which finally determines the actual 
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FIGURE 12 I Linear and quadratic responses of unit impulses (Ul) to 
pre-resonance stimulation frequencies {fi = 0.3, f2 = 1.3 and 
fs = 3.5 Hz). (A) Unit impulses corresponding to each action potential event 
of Figure 11. The dashed and dotted lines are the reconstructed linear and 
quadratic components based on the analysis of the second row. (B) Linear 
and quadratic analyses computed from the Fourier transform of the unit 
impulses. (C) The solid line is a slightly Gaussian filtered instantaneous 
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frequency (IF) trace constructed from the spike times. The dashed and solid 
lines are the reconstructed linear and quadratic components based on the 
analysis of the fourth row. The Gaussian filter is sufficient to remove peaks 
related to spike times. The analysis of less filtered IF curves (not shown), 
such that distorted spike events are visible, is very similar to scaled unit 
impulses analyses. (D) Linear and quadratic analyses computed from the 
Fourier transform of the filtered instantaneous frequency. 



unit impulses behavior, albeit in combination with the linear 
behavior. 

4. DISCUSSION 

The firing properties of the MEC stellate cells play a critical role 
in their fianction as part of the grid cell network. These neu- 
rons show oscillatory and nonlinear properties that are likely 
to be involved in the operation of networks involved in spa- 
tial awareness. We have used a novel theory, quadratic sinsoidal 
analysis (QSA), to rigorously determine the nonlinearity of SCs 
(stellate cells) near threshold for a direct comparison with lin- 
ear behavior. This multi-sinusoidal frequency probing not only 
provides a quantitative measurement of these properties through 
the QSA matrix, but also an algebraic characterization of the 
quadratic function as a Hermitian operator. The nature of the 



quadratic responses are significantly different than subthreshold 
linear behavior and gives an indication of the distinctive differ- 
ences between sub- and suprathreshold responses in SC neurons 
(Haas et al., 2007). We have shown that at physiological lev- 
els of stimulation, neurons and their models generate significant 
responses at harmonic and interactive frequencies that are not 
present in the input signal. Thus, the quadratic responses con- 
tain more frequencies over a wider frequency band than the input 
signal. As a consequence, they provide significant amplification at 
dynamically changing membrane potentials. 

Previous studies have shown that the nonlinear responses 
measured in neurons appears to be dominated by the dendrites 
(Magnani and Moore, 201 1), which receive the bulk of the synap- 
tic input. It is likely that the principal source of nonlinearity 
in stellate neurons near threshold is the persistent stochastic 



Frontiers in Cellular Neuroscience 



www.frontiersin.org 



August 2014 | Volume 8 | Article 239 | 14 



Magnani et al. 



Nonlinear analysis of entorhinal neurons 



Membrane potential -59mV (±3) 



Z (mV/nA) 




Ul Linear 




20 
18 
16 
14 
12 



R (mV/nA^) 



Eigenvalues (mV/nA^) 




0 10 20 

Unit Impulses (Ul) -56mV (i3) 



Ul R 



Ul Eigenvalues 




Hz 




10 



20 



-10 



QSA matrix Q 
(mV/nA^) 



r={-N,...,-i,+i,...+N} 



Ul QSA matrix 



I 



200 
150 
100 
50 




r={-N,...,-1,+1,...+N} 



FIGURE 13 I Linear and quadratic analyses of membrane potential and 
unit impulses for stimulation frequencies near the resonance {fi = 4, 
fz = 6.4, fs = 8.7 and fn = 12.1 Hz). (A) Analysis of subthreshold membrane 



potential. (B) Analysis of suprathreshold unit impulses. The linear and 
quadratic behaviors for the subthreshold and suprathreshold responses have 
similar form for these resonance frequencies. 



sodium conductance (^Nap)) which has been shown to be the 
primary source of channel noise in these neurons (Dorval and 
White, 2005). Thus, dendritic processing plays an important role 
in the firing behavior of stellate neurons in a number of ways 
(Zhuchkova et al., 2013). For example, the inherent filter prop- 
erties of cable structures in combination with their active voltage 
dependent channels alter the profile of the stimulation reaching 
the soma. We have shown that dendritic filtering greatly alters the 
nonlinear membrane potential response at the soma, which will 
depend on both the spatial location and frequency content of the 
stimulation. Finally, we have demonstrated that quadratic anal- 
ysis of modulated spiking behavior has significant similarities to 
membrane potential nonlinearities when the membrane potential 
is just below threshold. 

In conclusion, this analysis indicates that subthreshold lin- 
ear and nonlinear responses are similar to suprathreshold firing 
behavior. Thus, the combined linear and nonlinear behaviors 
near threshold of the membrane potential are reasonable esti- 
mates of suprathreshold behavior given by spike frequency mod- 
ulation. Linear and nonlinear behaviors a few millivolts below 
the threshold membrane potential are quite different with the 
nonlinear component being significantly reduced. The linear 
components of stellate neurons, such as resonance, are present 
at membrane potentials hyperpolarized to threshold despite the 
diminishing nonlinearities. In these membrane potential ranges, 
the linear resonance behavior is essentially due to the H current 
(Giocomo and Hasselmo, 2009). There is also a nonlinear compo- 
nent of the H conductance that can be measured at hyperpolar- 
ized values, however it is much smaller than nonlinear effects of 
the sodium conductances near threshold (Magnani et al., 2013). 



An additional effect of H conductances could occur if they are 
present in the dendritic tree. In this case, distal inputs could show 
bandpass resonance characteristics that would propagate to the 
soma involving both active and passive dendritic filtering, and 
further increase a nonlinear soma response. 
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